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ABSTRACT 

Development of new material models for describing the behavior of real materials 
(e.g. soils, ceramics and metals (monolithic and composite)) represents an important area 
of research in various engineering disciplines. This is evidenced by research activities 
associated, for example, with advanced high temperature metal matrix composites, 
reinforced concrete and geotechnical materials, to name a few. In particular, in 
applications involving elevated temperatures, the derivation of mathematical expressions 
(constitutive equations) describing the material behavior can be quite time consuming, 
involved and error-prone. Therefore intelligent application of symbolic systems to 
facilitate this tedious process can be of significant benefit. 

Presented here is a problem oriented, self contained symbolic expert system, named 
SDICE, which is capable of efficiently deriving potential based constitutive models in 
analytical form. This package, running under DOE MACSYMA, has the following 
features: i) partial differentiation (chain rule), ii) tensor computations (utilizing index 
notation) including both algebraic and calculus, iii) efficient solution f f r 

equations, iv) automatic expression substitution and simplification, v) back substitution of 
invariant and tensorial relations, vi) the ability to form the Jacobian and Hessian matrix 
and vii) a relational data base. Limited aspects of invariant theory ha.ve also been 
incorporated into SDICE due to the utilization of potentials as a starting point and the 
desire for these potentials to be frame invariant (objective). 

The uniqueness of SDICE resides in its ability to manipulate expressions in a 
general yet pre-defined order and simplify expressions so as to limit expression g rc >wt . 
Results are displayed, when applicable, utilizing index notation. SDICE has been desi 8" ed 
to aid and complement the human constitutive model developer. A number of examples 
will be utilized to illustrate the various features contained within SDICE. It is expected 
that this symbolic package can and will provide a significant incentive to the development 
of new constitutive theories. 


1. INTRODUCTION 

Efforts in constitutive research involve for example, the development of 
mathematical relationships for predicting reversible and irreversible material response, 
derivation of material stiffness matrices appropriate for finite .element calculations, 
characterization of model parameters, formulation of a Jacobian matrix for implicit 
numerical integration schemes, issues of convexity, and computer implementation, me 
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entire process requires significant manual algebraic manipulations and computer 
programming. Hence, the response time for the related efforts is quite long. As a result, it 
can be rather difficult to introduce significant changes or modifications into a constitutive 
theory. Moreover, the outcome of the research effort may be affected by human errors 
which are often difficult to detect. In this regard, symbolic computation can play a major 
role. Furthermore, application of symbolic manipulation can provide a significant 
incentive to the development of new constitutive theories and their applications, e.g. finite 
element. Presented here is a description of a software package named SDICE (Symbolic 
Derivation of Constitutive Equations) which was designed to aid and complement the 
human constitutive model developer in constructing, analyzing and implementing potential 
based constitutive equations. 

Results obtained from direct application of a general purpose symbolic system, such 
as MACSYMA (see Ref. 7), have been shown to be not useful in most cases due to the 
number of steps involved in the derivation process and the problems associated with 
expression growth (Refs. 4,5 and 17 — 22). For this reason, resourceful derivation procedures 
(which are continuously in communication with a relational data base) must be developed 
so that optimal results can be obtained. This is the intent behind the development of 
SDICE, whose uniqueness resides in its ability to manipulate expressions in a pre— defined 
order and simplify expressions so as to limit expression growth. The essential features of 
the approach taken in SDICE to address the above problems consists of: 

1) a structured derivation procedure to avoid redundant steps and to minimize 
expression growth, 

2) implementation of special procedures (e.g. procedures for simplification and 
pattern match) to facilitate the derivation process, 

3) expression substitution and simplification during the entire derivation 
process by incorporating several levels of processing, 

4) automatic grouping and labeling of common factors (intermediate variable 
substitution), 

5) taking advantage of permutation and symmetry relationships of terms 
involved in each derivation step, and 

6) as a rule— based system, intended for constitutive equation research, SDICE 
will record user defined rules and store them in a relational data base, 
whereby the information may be retrieved, redefined and restored as 
required. 

These six features are interconnected in such a way as to allow them to work in concert 
with one another, thereby providing the user with simplified final expressions. In the next 
section a brief description of potential based constitutive equations is given. This section is 
followed by an overview of the special problem oriented functions presently available in 
SDICE as well as a section illustrating their application to a transversely isotropic 
formulation. 


2. APPLICATION TO POTENTIAL BASED CONSTITUTIVE EQUATIONS 

Constitutive laws provide the link between stress components a - and strain 

components e^ at any point in a body. These laws may be simple or extremely complex, 

depending upon the material of the body and the conditions to which it has been subjected. 
Consider the well known case of a hyperelastic material. Here, the stress and strain 
components are related through a normality structure utilizing either a strain energy or 
complementary energy function, i.e. 


or 
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e - dQ 
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For inelastic material behavior the internal state variable potential viewpoint is 
adopted, i.e., 


ft — ft((7jj,a^,T) 

with the generalized normality structure 


( 3 ) 


and 
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<* 0 = -h(a jan 


0,7= 1,2,... ,n 


(5) 


as discussed in Refs. 10 through 12. Where Q is the complementary dissipation potential 
function, e.j the inelastic strain, a q the internal state variables, n the number of internal 

state variables and (') denotes differentiation with respect to time. Equations (4) and (5) 

are known as the flow and evolutionary equations, respectively. 

Frame invariance (objectivity) of the resulting constitutive relations is insured by 
requiring the potential, ft or W, to depend only on certain invariants and invariant 
products of its respective tensorial arguments, i.e., an integrity basis, see Ref. 16. Both 
isotropic and transversely isotropic material symmetries have been considered thus far. 
Transversely isotropic material symmetry is included in the potentials of Eqs. (1) — C3) by 
introducing a directional tensor u-m, e.g. ft=ft(<j.j,a^,u.Uj,T) or W=W(e-j,u-Uj). The 

symmetric tensor u^Uj is formed by a self product of the unit vector u. which denotes the 

local preferred direction. 

Two viscoplastic theories have been proposed for isotropic materials, see Refs. 13 
and 14, and employed to verify implementation of various special functions within SDICE 
as described in Refs. 1,2 and 23. In both theories the existence of a dissipation potential is 
assumed, and the form is taken to be, 




( 6 ) 


where the dependence of the applied stress and internal stress (cf. Eq (3)) enters through 
the scalar functions F^j) and G(a i j), respectively. 

For a material with transversely isotropic symmetry the dissipation potential is 
assumed to take the form of Eq. (6) where the dependence of the applied stress, internal 
stress and preferred direction enters through the scalar functions F(S-j,UjUj) and G(a|j,u-Uj) 

respectively. The stress dependence is given by 


where 


= [AJ 2 + BJ 5 + CJ 4 ] - 1 

(7) 

= [AJ 2 + BJ 5 + cj2] 
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ij ij j kk ij 


Upon application of Eqs. (4) and (5) the resulting flow and evolutionary laws are: 

r. 


i = f(F) 


»ij = h (G)'y - -KGJify 

respectively, where 

r ij = AE ij + B I u k u i S jk + u j u k E ki -^- J 4 u i u j] + ^ CJ 4< 3u i u j - V 


( 10 ) 

( 11 ) 


"ij = Aa rj + B t u k u i a jk + u j u k a ki - y J 4 u i u j) + i CJ 4 ( 3u i u j - J ij) 

The above transversely isotropic equations will be employed in subsequent examples to 
illustrate specific functions incorporated in SDICE. 

3. SPECIAL PROBLEM ORIENTED FUNCTIONS 


In deriving material constitutive equations and matrices, a number of general 
mathematical capabilities are desirable, including: 

1) partial differentiation (of both scalar and tensorial functions), 

2) tensor computations (utilizing index notation), 

3) matrix solution of systems, 

4) factorization of common terms, 

5) expression substitution and simplification (i.e. intermediate scalar and 
tensorial variables), 

6) back substitution of invariant and tensorial relations, 
as well as some specialize capabilities such as 

7) formation of the Jacobian and Hessian matrix 

8) solution of eigenvalues, 

9) the ability to input and apply simplifying conditions, 

10) the ability to convert from index notation to matrix notation, 

11) the automatic determination of an integrity basis (invariant theory), 


4 



12) the automatic fortran generation of the resulting expressions, 

and last, but not least, the desire to perform these various tasks within a user friendly 
environment. 

Symbolic packages, such as MACSYMA, MAPLE, REDUCE, MATHEMATICA, 
etc., have varying degrees of capability to address certain aspects of the above 
requirements, however, direct use of a symbolic system in the derivation process can be 
complicated and lead to rapid expression growth. Therefore, special purpose procedures 
have been and are being designed and implemented into SDICE in order to simplify and 
expedite the derivation process. These special procedures are implemented as separate 
LISP (Ref. 24) functions, thus enabling manipulations of the internal data structures of the 
MACSYMA expressions as well as compensation for the inefficiency and lack of 
simplification mechanisms within the existing MACSYMA level built-in functions. 

For example our work thus far in dealing with the differentiation and the 
minimization of expression growth of implicit tensorial functions has resulted in the 
development of several strategies, namely: 

1) Store the invariant and tensorial relations in a relational data base. 

2) Implement procedures to compute tensor expressions according to the rules 
defined in tensor calculus. 

3) Map a tensor into the domain of a scalar by utilizing a property list to store 
the variable and its subscript; thereby, allowing all differentiation to be 
treated in the same way. 

4) Generate subscripts for intermediate tensor variables and store them in the 
same property list in a pre— defined order. 

5) Represent differentiation results by a search tree, starting with the potential 
function as the root of the tree and its descendants without dependent 
relations among themselves as leaves. A separate procedure decides whether 
the function and its variables are tensorial or scalar. 

6) Finally, check the property list and if the function involves tensors, 
subscripts are added back for the final result according to the predetermined 
order. 

7) Simplification of the result is accomplished by i) checking the data base so as 
to replace any known invariant or intermediate relation with its 
corresponding name in the final expressions, ii) grouping common terms by 
factorization and iii) identifying and naming terms that can be written as a 
tensorial variable. 

These strategies as well as others, are contained within a new set of problem 
oriented functions residing in SDICE which supplement those normally available under 
MACSYMA. It is important to note that these functions utilize some relevant 
MACSYMA functions, such as EV, RATSIMP, FACTOR, etc. However they are not 
merely driver routines, but functions with problem oriented computational algorithms 
embedded in them. A list of the name, argument list and short description ofthese new 
problem oriented functions is given in the Appendix. 

4. APPLICATION OF SDICE FUNCTIONS 

Here a number of the available functions within SDICE will be exercised. Input 
required by SDICE is entered interactively through a user friendly interface and stored 
within a relational data base. SDICE provides the user with the ability to create, delete, 
modify or project a given data base as well as add, check or delete common rules, through a 
data base manager. Due to the fact that MACSYMA typically returns lower case letters 
and the desire to utilize the greek alphabet, the following notation has been adopted within 
SDICE: 
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a — z 
aa — zz 
ga-gz 
gaa— gzz 


lower case English alphabet 
upper case English alphabet 
lower case Greek alphabet 
upper case Greek alphabet 


Also note that commas between indices are merely separators and are not to be associated 
with spatial derivatives as standard index notation suggests. 


4.1 Invariant Theory 

Utilization of potential functions as a starting point and the desire for these 
potentials to be frame invariant (objective) prompted the introduction into SDICE of 
limited aspects of algebraic invariant theory. As a result, given the number of symmetric 
second rank tensors contained within a potential and any user defined simplifying 
conditions, SDICE will determine automatically the corresponding integrity basis. For 
example, consider a function which contains two second rank symmetric tensors, defined as 


and 


U U U 


and one simplifying condition u.u- = 1. If we invoke SDICE by employing the function 
called INV, 


in20 ==> inv ( [gss, uu] , [uu*u [i ] *u ( j 3 / u [i ] *u [i ] *1 ] ) ; 

the following result is obtained. 

out 20 ==> [[[css], gss ], [[gss ] , gss gss ], 

i, i i# j 3/ 1 

3 - 

f [gss ], gss gss gss j, [[gss uu, gss uu ] , u gss 

- » j j » k k, i i 

2 2 2 

[ [gss uu, gss uu ] , u gss gss u ] ] 

i i, j j, k- k. 


Theoretically (see Ref. 16) the integrity basis for a function containing two second rank 
symmetric tensors should be comprised of ten invariants. However due to the special 
nature of uu, three invariants are found to be equal to one and are therefore discarded, 
while two other pairs of invariants are found to be equivalent to one another. Thus 
reducing the remaining seven by two. Therefore due to the imposed condition, u^Uj = 1, 

the above five invariants are seen to comprise an integrity basis for this potential with 
these two symmetric second rank tensors. Note that brackets enclose a specific invariant 
denoted using index notation while sub— brackets define the invariant using matrix 
notation; for example the first invariant is interpreted as 

[gss] = trE 


in matrix notation or 


6SS i,j $ % 


in index notation, while the fourth is 
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trEu = trEu z or 


u.E. a* • 

i U J 


4.2 Partial Differentiation 

Determination of flow and evolutionary laws are directly linked to the ability to 
perform partial differentiation as discussed above and illustrated in Eqs. (4) and (5). This 
capability is accessed within SDICE through the function named DIF (see the appendix for 
a description of the argument list). For example considering Eqs. (4) and (6) through (9), 
it is apparent that in deriving the flow law for the above transversely isotropic model, the 
inelastic strain rate is given by: 


iJ 


= du 



dSl 


'dF dJ 2 + dF 3J 5 +dF dJ A ' 

aE kl SS r s 

<?F 


SJ 2 5J 5 3J 4 cTEy 



Notice that three types of partial derivatives must be taken, i) one in which both the 
function (numerator) and variable (denominator) are scalars, ii) one in which both the 
function and variable are tensors and iii) one in which the function is a scalar and the 
variable is a tensor. All three types are included in the DIF function. Upon employing the 
data base manager to store Eqs. (6) thru (9) in a data base named anij2, the above flow 
law may be obtained by issuing tne following command. 


ini «> dif (goo, gs [i, j ] , ani j 2 ) ; 


The result returned is then: 

2 

out l ==> gk <f(ff) (6 cc u u + (- 2 cc - 2 bb) gd ) jj4 

i j i, j 

+ f(ff) (u (3 bb u gss + 3 bb gss u ) + 3 aa gss ))/(6 gm) 

il i j, i 1 i / il j ^ 3 

Intermediate tensorial variables can then be established by employing the function, 
GROUPT, as shown below; 


in2 « > ed [i, j ] : groupt (% , [i, j ] ) ; 

2 

out2 ==> gk: (f(ff) (6 cc u u + (-* 2 cc - 2 bb) gd ) j j 4 

i j i# j 

+ f(ff) (3 bb ggtl + 3 aa gss )>/(6 gm) 

i > j i' 3 


in3 ==> ggtl [i , j ] ; 

out3 ==> u u gss + gss u u 

i il j, il i, il il j 
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Here only one level of tensor substitution (intermediate variable substitution) is invoked. 
Comparing the expression denoted by out2 above with Eq. (10) we see, after minor 
simplification, that they are equivalent. This final factoring out of common terms such as 
f(ff) and finding the minimum common denominator is still a topic of ongoing research. 
Additional details regarding the differentiation algorithm can be found in Reis. 2 and 23. 

4.3 Eifienvalue Problem 

A question of importance to constitutive developers is whether or not the flow (or 
yield) surface one is constructing or employing is convex. One way to ascertain this is to 
calculate the corresponding eigenvalues and examine their sign. Here we will determine 
the eigenvalues corresponding to the function F, see Eq. (7). This example was selected 
based on the fact that numerous SDICE functions are utilized, e.g. INPCON, GETF, 
APPCON, RATSMP_A, FORMM and CHAR. 

Let us begin by imposing specific conditions on the function F, i.e. assume a 
principle stress state, take the internal stress state to be zero, and examine a specific 
direction orientation, u*. This is easily accomplished by invoking the function INPCON. 

in5 «**> inpcon(l); 


TEST 1 

PLEASE INPUT THE CONDITIONS (type d or done when finished) 
ENTER THE CONDITION : gs [ i , j ] «= C tP / 0. , 0 1 , [ 0 , ? , 0 ’ [ 0 , 0 , p ] ) ; 
ENTER THE CONDITION : ga [ i , j ] = [ [ 0 , 0 , 0 ] , [0,0,0, [0,0,0]]; 
ENTER THE CONDITION : u [i ] * [ 1 , 0 , 0 ] ; 

ENTER THE CONDITION : d; 
out5 ==> done 


Equation (7) is retrieved from the data base, anij2, by invoking the GETF function. 


in6 =*=> getf(ff,anij2); 

2 

out6 ==> bb j j 5 + cc j j 4 + aa j j 2 - 1 

The set of conditions previously entered with the INPCON command are then applied to 
the above expression using the APPCON function. 

in7 ==> appcon ( [ 1 1 , % ) ; 

2 2 
eqnl : { (gs + (2 gs - 4 gs ) gs + gs - 4 gs gs 

3 , 3 2 , 2 1, 1 3, 3 2 , 2 1, 1 2, 2 

2 2 2 
+ 4 gs ) cc + (gs + {2‘gs - 4 gs ) gs + gs 

1,1 3/ 3 2, 2 1, 1 3, 3 2, 2 

2 2 

- 4 gs gs + 4 gs ) bb + (3 gs + (- 3 gs - 3 gs ) gs 

1, 1 2, 2 1, 1 3, 3 2, 2 1, 1 3, 3 

2 2 
+ 3 gs - 3 gs gs + 3 gs ) aa - 9) /9 

2 , 2 1 , 1 2 , 2 1 , 1 
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out 7 -*> [eqnl ] 

Due to the scalar nature of F, see result out6, only a single equation is returned. If the 
expression of interest were tensorial then multiple equations could be returned depending 
upon the imposed conditions. The diagonal components an< ^ ^33 ^ ^ en 

replaced by their principal stress representations and a ^ respectively, by employing 

the alternative form of the APPCON function. 


in 8 “> appcon (eqnl, [gs [1 , 1 ] =gsl, gs [2, 2] =gs 2 , gs [ 3 , 3 ] «gs 3 ] ) ; 

2 

out8 ==> ( (cc + bb + 3 aa) gs 3 + ({2 cc + 2 bb - 3 aa) gs 2 

2 

+ {- 4 cc - 4 bb - 3 aa) gsl) gs 3 + (cc bb + 3 aa) gs 2 

•r (- 4 cc - 4 bb ’ 3 aa) gsl gs 2 + (4 cc + 4 bb + 3 aa) gsl - 9 ) / 9 


The nonphysical coefficients A,B, and C are now replaced by their physically 
meaningful counterparts, , 7 / and u (denoted by gk, ge, go, respectively) determined 

elsewhere, by again utilizing the function APPCON. 

in9 ==> appcon { % , [aa*l/gk."'2, bb=* - (ge~ 2 - 1 ) / (ge A 2*gk.~2) , 

cc sa l/gk''2 A ( ( 2 + go ~ 2 ) / ( 4 * gc A 2 - 1 ) - I . ae' 2) ] ) ; 


2 2 2 2 2 2 2 
out9 ==> (- 4 gk go + gs 3 (gs 2 (1 - 2 go ) - gsl) + go gs 3 + go gs2 

2 2 2 2 2 

- gsl gs2 + gsl + gk ) i i 4 gk go - ck ) 


The resulting expression is then simplified utilizing the RATSMP__A procedure. This 
function differs from MACSYMA’s simplification function RATSIMP in two primary ways. 
First, it allows a priority by variable or expression to be specified during the simplification 
process, thereby giving the user control over the appearance of the final expression. 
Secondly it performs intermediate variable substitution operations throughout the 
simplification process, thus resulting in more compact and tractable final expressions. 


inlO ==> ratsmp_a ( 5, [gsl , gs 2 , gs 3 ] / 1 ) ; 


qql 


2 

go 


2 

gk (2 go - 1) (2 go + 1) 


2 

2 go - 1 

cq2 = 

2 

gk (2 go - 1) (2 go + 1) 
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1 

qq3 

2 

gk (2 go - 1) (2 go + 1) 

2 2 2 

outlO «“> qql gs3 + qql gs2 + qq3 gsl - gsl gs3 qq3 - gsl gs2 qq3 

- gs2 gs3 qq2 - 1 

At this stage one can either employ the function EIGEN, described in the appendix, 
or make use of the various procedures utilized by EIGEN to perform the task. In order to 
give the reader a clearer understanding of the various features of SDICE, the latter option 
has been selected. 

The following quadratic polynomial can be written in matrix form by utilizing the 
FORMM function, that is 

inll "> f ormm { % , [gs 1 , gs2 , gs3 1 ) ; 
out 11 “*> x aa x - 1 

i i, j j 

where 

i n 1 2 ==> x [ i ] ; 

outl2 *=-> [ gsl gs2 gs3 ] 


in 1 3 ==> aa (i, j ) ; 

C qq3 qq3 ] 

C qq3 ] 

[ 2 2 ] 

C ] 

[ qq3 qq2 ] 

out 1 3 = = > [ qql ] 

[2 2 ] 
[ ] 

( qq3 qq2 ] 

C qql ] 

[22 ] 


i n 1 4 = = > x [ j ] ; 

[ gsl ] 
[ 1 

out!4 [ gs2 ] 

[ ] 

[ gs2 1 


The characteristic equation is then determined using the CHAR procedure. 


in 15 =-> char {a a [ i , j } ) ; 


out 15 ==> - (-2 qqz + qq2 + 2 qql' i. - 2 cgz r qq3 + 2 ggz qq3 + qq2 qq3 

- 2 qql qq3 - ggz qq2 + 2 ggz qql ) / 4 
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where ggz represents the typical eigenvalue A. The characteristic polynomial is then 
simplified using intermediate expressions as follows: 


ini 6 ==> ratsmp_a ( % , ggz, 1 ) ; 

qq4 = (qq2 + 2 qql) qq3 (qq3 + qq2 - 2 qql ) 

2 2 2 
qq5 = 2 qq3 - 8 qql qq3 + qq2 - 4 qql 


qq6 = qq3 + 2 qql 


3 2 ggz qq5 qq4 

out 1 6 ==> - ggz + ggz qq6 + 

4 4 


inl7 — > appcon(qq4); 
out 1 7 -=> 0 

inl8 ==> outl6+' qq4 /4; 

3 2 ggz qq5 

out 1 8 = = > - ggz + ggz qq6 + 

4 


The eigenvalues are finally determined by using the MACSYMA SOLVE function, i.e. 

inl9 ==> solve (%, ggz); 


outl9 **=> [ggz 


2 

qq6 - sqrt (qq6 + qq5) 


2 


2 

sqrt Cqq6 + qq5) + qq6 
ggz , ggz = 


Inequalities can then be generated, if necessary, to restrict the material parameters so that 
the surface is convex. Development of procedures to assist in the establishment of these 
inequalities is an extremely challenging task and is an area of current research. 

5. CONCLUSION 

A problem oriented, self contained symbolic expert system, named SDICE, which 
runs under DOE MACSYMA and is capable of efficiently deriving a special class of 
constitutive equations in analytical form, namely those derivable from a potential function, 
has been presented. A new set of special purpose procedures have been developed which 
provide capabilities in areas such as; partial differentiation (chain rule) of both scalars and 
tensorial functions, efficient solution of sparse systems of equations, automatic expression 
substitution and simplification, back substitution of invariant and tensorial relations, 
formation of Jacobian and Hessian matrices, and algebraic invariant theory. A number of 
examples have been included to illustrate various features contained within SDICE. 

The uniqueness of SDICE resides in its ability to simulate the human intelligence 
required to manipulate expressions in a general yet pre— defined order and simplify 
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expressions so as to limit expression growth. This simulation is accomplished though the 
utilization of a relational data base, in that all problem— oriented algorithms are 
continuously in communication with this relational data base thus allowing automatic 
expression and variable substitution throughout the solution process. 

Work is continuing, with emphasis upon 1) further improvement in the 
simplification capabilities available, 2) in a user friendly interactive interface, 3) automatic 
Fortran code generation of the resulting expressions and 4) automation of the identification 
of inequalities. In summary the application of intelligent symbolic methods in the 
development, analysis and implementation of constitutive theories has thus far shown 
widespread benefits in expediting and reducing the error— prone nature of this tedious 
process. Public release of SDICE is anticipated in the near future. 
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APPENDIX - LIST OF AVAILABLE SDICE FUNCTIONS 


APPCON(fnl,n2,..],obj) 

Applies [nl,n2,...] test conditions to the object obj, where obj can be either an atom 
or a tensor with appropriate indicies. 

Alternate form 

APPCON([eqnl,eqn2,..], [conditions]) 

Simplifies eqnl,eqn2,etc. corresponding to the specified conditions given in the 
second argument list. 

CHAR(m) 

Calculates the characteristic polynomial of the matrix m. Intermediate variables 
are used to simplify the results. 

DIF(fnl,fn2,dbn) 

This function differentiates fnl with respect to fn2 , where fnl or fn2 may be either 
a scalar or tensor expression. If fnl is an implicit function, DIF will automatically 
apply the chain rule to the extent required. Similarly, the form of the result is 
consistent with the initial form of fnl and fn2 defined in the data base. Both fnl 
and fn2 are defined in the data base dbn. 

EIGEN(exp,varl,var2,..) 

Calculates the eigenvalues of the matrix formulated from the given quadratic 
expression ,exp, with dependencies varl,var2,etc. Combines functions FORMM and 
CHAR into one command. 

FORMM(exp,[varl,var2,..]) 

Formulates the matrix form of the expression, exp. Intermediate variables are used 
to simplify the results. 

GET_LEAD_M(m,s) 

Obtains the leading principle matrix m with dimensions s by s. 

GETF(fn,dbn) 

Obtains the definition of the function fn defined in the data base dbn. 
GROUPT(exp,(k,l]) 

Groups tensor product terms, possessing common subscripts [k,l] and automatically 
generates a tensor intermediate variable ggti^ j which represents these terms. 

HESS([fnl,fn2,...],[varl,var2 J ...],dbn) 

Calculates the Hessian matrix. Note fnl,fn2,etc. are defined in the data base dbn 
with varl,var2,etc. as their dependents. 

INPCON(n) 

Allows user to enter n sets of test conditions, which will be used by the APPCON 
function. These conditions remain available until one exits from the execution 
sub— level of SDICE. 

IN V([tenl,ten2,...], [conditions...]) 

Automatically derives the integrity basis for a given function, provided the tensors 
(tenl,etc.) are symmetric and of rank two. The second argument can contain the 
definition of the tensors as well as any other condition or simplification rules one 
wants applied. 

JAC([fnl ) fn2,...],[varl,var2,...],dbn) 

Calculates the Jacobian matrix. Note fnl,fn2,etc. are defined in the data base dbn 
with varl,var2,etc. as their dependents. 
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LINSOLVE2([eqnl,eqn2,...J,[varl,var2 ) ...]) 

Solves the system of linear equations simultaneously for the list of variables. This 
function is significantly more efficient than the existing package under MACSYMA 
in handling sparse systems. 


PUTF(exp.dbn) , _ . . 

Stores the expression exp as the definition of a function (named by the user) in the 

data base dbn. 


RATSMP(exp vl J v2 ) .. J vn)^ : 

Enables rational simplification of the expression exp relative to the specification of 
the variable (function or expression) vn. The order is important since vl has the 
highest priority and vn the lowest. If a variable within exp is missing then this 
variable is given lower priority then the rightmost vn. 

RATSMP_A(exp,vl > v2,..,vn,flag) : * : 

Enables rational simplification of exp relative to the priority of vl,v2,...,vn, 
respectively. However with this function call intermediate variables are generated 
to represent sub expressions of variables of lower priority and are displayed if flag is 

* Note in both RATSMP functions multiple variables can be given the same priority by 
enclosing within brackets. 

SETDBN(dbn) 

Identifies dbn as the working data base. 

SETPRI(varl,var2,...varn) . . 

Defines the global priority of variables, with respect to simplification routines. Varl 
has the highest priority while the others follow sequentially. 


SHOWPRI 

Shows the current priority list. 


SOLVE2([eqnl,eqn2,...],[varl,var2,...]) < 

Solves the list of simultaneous polynomial equations (linear or nonlinear) for the list 
of variables. It utilizes the SOLVE function provided by MACSYMA as well as the 
simplification functions provided in SDICE. 


VWCON(n) 

Displays the nth set of test conditions. 
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